Transitions at avoided level crossing with interaction and disorder 
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We investigate the influence of interaction between tunneling particles and disorder on their 
avoided-level-crossing transitions in the fast-sweep limit. Whereas the results confirm expec- 
tations based on the mean-field arguments that ferromagnetic/antiferromagnetic couplings sup- 
press/enhance transitions, we found large deviations from the mean-field behavior for dipole-dipole 
interactions (DDI) in molecular magnets Mni2 and Fes. For ideal crystals of the needle, spherical, 
and disc shapes DDI tends to enhance transitions. This tendency is inverted for the needle shape 
in the presence of even small disorder in the resonance fields of individual particles, however. 
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PACS numbers: 03.65.-w, 75.10.Jm 

Transitions at avoided level crossing or the Landau- 
Zener (LZ) effect |l|,|2j is a well known quantum phenom- 
ena, mainly in physics of atomic and molecular collisions. 
In the time-dependent formulation the dynamics is that 
of a two- level system described by a pseudospin s = 1/2 
whose Schrodinger equation (SE) can be solved exactly 
for the linear sweep of the effective longitudinal mag- 
netic field |2|. As the SE for a spin 1/2 is mathematically 
equivalent to the classical dissipationless Landau-Lifshitz 
equation (LLE), the effect can be envisioned classically 
as a rotation of a magnetization vector. 

Recently the LZ effect was observed in the solid-state 
world on crystals of molecular magnets Fes (see Refs. 
0,01 for a recent review). This posed a new major prob- 
lem of the LZ effect in many-body systems with interac- 
tions. As tunneling of one particle between the two states 
changes conditions for the others [in this case mainly 
via dipole-dipole interactions (DDI)] one is confronted, 
in general, with a SE for N 1 coupled two-level sys- 
tems that contains 2 N time-dependent coefficients. 

The problem can be simplified if one applies the mean- 
field approximation (MFA) that considers one particle 
tunneling in the effective field that is a sum of the ex- 
ternally sweeped field and the molecular field from other 
particles that is determined self consistently. This is a 
model of the nonlinear LZ effect that was applied to tun- 
neling of the Bose-Einstein condensate 0, 0] ■ Again this 
problem can be reformulated in terms of a classical non- 
linear LLE. The MFA solution shows that ferromagnetic 
interactions suppress transitions while antiferromagnetic 
interactions enhance them. Corrections to the MFA for a 
simplified "spin-bag" interaction model of TV pseudospins 
coupled to all others with the same coupling J 0,0 were 
studied in Ref. jldj . However there are no rigorous re- 
sults for competing interactions such as the DDI, while 
existing theories use postulated rate equations (see, e.g., 
Ref. 0). 

The aim of this Letter is to develop, for arbitrary in- 
teractions, an expansion of P for fast sweep rates v where 
P is close to 1 . Whereas the first term ~ 1 /?; in 1 — Pis 
in most cases insensitive to the interaction (that is why 
it was possible to extract the correct value of the ground- 
state tunnel splitting A in Fes 0)) the next term ~ 1/v 2 



does depend on the interaction and it shows whether 
transitions are enhanced or suppressed. We take into 
account inhomogeneities of individual resonances as they 
can strongly reduce the effect of interaction. Then we ap- 
ply our result to the DDI and show that its effect differs 
considerably from the mean-field prediction. 
We consider the transverse-field Ising model 
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{[H z (t)-Vi\ 



. (1) 

where <Tj are Pauli matrices, H z {t) — vt is the time-linear 
sweep field, Vt is the local shift of the resonance field, A 
is the splitting of adiabatic energy levels for = 0. The 
initial state of our model is all pseudospins down. For 
Jij = the well known solution [lj, |2| for the final-state 
probability for a spin to remain in the initial state is 



P = P(oo) = e" 



ttA 2 
2^7' 



(2) 



whereas P(t) can be expressed via hypergeometric func- 
tions. The fast-sweep expansion of Eq. @ is P = 



1 



2 /2- 



We will see that interaction J,-, mod- 



ifies the term e 2 . At order e 2 it is sufficient to take into 
account maximally two spin flips out of the initial state 
and to write the wave function in the form 

*(*) - C (t) III ... J,) + c i( t ) (7 i+ |U • • • I) 



h^ J Cij{t)vi+V3+ III ••• I) 



(3) 



The one-particle staying probablity averaged over the 
sample is 



(4) 



The Schrodinger equation reads 
A 



ihdo = x cq 



ihdt = Ei{t)ci - —cq - -^i 
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ihdij = E i3 (t)cij - — (c 



A^ 



(5) 



etc. Here E are the eigenvalues of the Hamiltonian with 
A = and the ground-state energy E (t) subtracted and 

Ei{t) = -H z {t)+V h V i = V i + 2^Jij 

3 

Etjit) = ^H^ + Vi + Vj-Uij. (6) 

For the fast linear sweep the solution of Eqs. 10 is a 
series in the integer and half-integer powers of e that can 
be solved by iterations starting from co(t) = 1. It is suf- 
ficient to retain cq, Ci, and while Ciji can be dropped. 
Inserting found coefficients into Eq. (@J and calculating 
double and triple time integrals yields the final result 



a reasonable precision within a reasonable time. In the 
homogeneous case 7^ = and 1^ simplifies to 

= c (/y [1 - c (/y ] + s (/y [1-5 (/yji3) 

The limiting forms of F{j3^j) are 



F(P«) = { Pi 



f*J cos (!/% + !), 

PiJ, |/?yl<l 
ft, » I- 



For the weak interaction Eq. Q) then yields 
P^ 1 -e+ - 



4Jn 



^A ' 



(14) 



(15) 



P^ 1 -e + e 2 /2 + e 2 I Q , 



(7) 



where 



^ + cos (2vr 7 g ) /3 lJ ) ^ + sin faff Pij) 



Zj,- = Iji, and Aj,-, fl,, Ch are defined by 
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\ ~ S (7ji) 
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2 - 5 



2 ~ 5 UrJ 



(8) 



(9) 



(10) 



a generalization of Eq. (26) of Ref. {9j for the arbitrary 
form of Jtj. Note that Eq. (|15|) is essentially a MFA re- 
sult as it only depends on the zero Fourier component Jo 
of the coupling . In contrast to thermodynamic sys- 
tems, here the applicability of the MFA is not controlled 
by the interaction radius alone. For the nearest-neighbor 
interaction with z neighbors the relative correction to 
the last term of Eq. ifTB) is - [4J /(ttA)] e 1 ' 2 /z. Eq. 
shows that ferromagnetic interactions, Jy > 0, increase 
P and thus suppress transitions, whereas antiferromag- 
netic interactions facilitate transitions. The saturation 
for strong ferro- and antifcrromagnetic interactions in Eq. 
(fH|l corresponds to the case of well-separated resonances 
studied in Sec. Ill of Ref. 9]. 

Let us proceed to the inhomogeneous case, 7^ ^ 0. 
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max(|/3 y |,l) > 1, Eq 

r ~ /"IT 2 7T 

7i . = _ cos ^_ 7 .. + _ 



yields 
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Here C{x) and S{x) are Fresnel integrals and 

Pij > 



7« 
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Eqs. lj7|- l|12H is our main result that is valid for arbi- 
trary interactions and resonance shifts. Note that it has 
a pair structure and thus it can be verified against the 
direct numerical solution for the model of two coupled 
particles. Analytical form makes its application practi- 
cally possible; Triple time integrals that arise at the in- 
termediate stage cannot be computed numerically with 



i.e., strong inhomogeneities reduce the effect of interac- 
tions as individual resonances are well separated and flip 
of one particle does not bring another particle past or 
before the resonance by the changing effective field. For 



(11) 1 0ij\ < 1 Eq. © yields 



. /7T (0)2\ . /7T 

(■rSO-tfOC) 



,(0)2 



Pi 



S 



(17) 



For 
i.e., 



7 (°) 
1 1] 



<C 1 this simplifies to 1^ = ft^ ^1 — ttj^ 2 /2 



weak inhomogeneities do not essentially suppress 



weak interactions. For |/3 i 



(0) 



3> 1 one obtains 




Pij < 
Pa > 0. 



(18) 



As the first line of Eq. I|18fl corresponding to strong 
antiferromagnetic coupling oscillates fast with 7^ the 



3 



limiting value F = —3/2 in Eq. (|14|) is unstable with 
respect to small inhomogeneities V*. The effect of inho- 
mogeneities can be accounted for by averaging Eq. © 
over stochastic values of on with a normalized Gaussian 
distribution p a (a) = (2tt(5 2 )~ 1/2 exp [-a 2 / (2S 2 a )] and 
quadratic average (a?) = 5 2 a . The distribution of 7^ = 

a% — ay is then given by the same function with 8 a 28 2 a . 
For 8 a <C 1 one can just set 7^ , 7 =>■ /3y in Ay , Bij , and 

Cy in Eq. © and use <Jcos (27r 7 g ) /3 u ) ) = e-( 27r/3 « 5 ") 2 . 
As this factor decays at /3y ~ 1/£ Q ^> 1, one can further 
simplify the result and replace of F{j3 io ) of Eq. 113fl by 
T(j3 i3 ) that satisfies F(±oo) = ±1/2: 



F(J3 iit 6 a ) = FtfiJ + SFfaM, 



8 a < 1 



_ P -i 2 -p^T 



For (3^ <C 1 from Eq. (|17|) one obtains 
F((3 l3 ,S a )= (3^/(2^1), f{x) 
x + 1 x — 1 



(19) 



1 vvi + ^-i 

Function f(x) monotonically decreases and satisfies 
1 - x/2, x<l 



(20) 



/(*) = 



1/V2x, x>l. 



(21) 



Let us now turn to the DDI between tunneling spins 
±5 of magnetic molecules aligned along the z axis: 



Jij — 



3 cos 2 % - 1 
«o 3^ , (22) 
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FIG. 1: G of Eq. 1)240 for the sphere vs width of distribution 
of individual resonances 8 a . Dashed lines on the left and right 
are asymptotes of Eqs. I|26[l and l|27|l . respectively. 



the result independent of the lattice structure. Since at 
large distances 7y ~ /3y behaves as the DDI, the result 
depends on the sample shape. For the model with ran- 
dom resonance shifts using Jy = i 7 '(/3 ij -, 5 a ) one obtains 



4tt/ (2k51) 



J = Gf, G = G (Sphcrc) + (l/3 

(24) 

where rS z ^ — 1/3, 0, and 1 for a sphere, needle and disc, 
respectively, f(x) is that of Eq. J2DJ, and 



^(Sphere) = Rf ^2 



K 




9^ 

1 ln \ / 3 + 1 ^ 
\/3 n \/3 - 1 1 



d/3 



F(0,6 a ) 
P 2 



-0.66924. (25) 



where vq is the unit-cell volume, ry is the distance be- 
tween the sites i and j, and cos#y = ry^/ry. We con- 
sider 5* = 10 molecular magnets Mni2 having a tetrago- 
nal lattice with parameters a = b = 17.319 A, c = 12.388 
A (c is the easy axis) and vq = abc — 3716 A 3 and Fes 
having a triclinic lattice with a = 10.52 A (a is the easy 
axis), b = 14.05 A, c = 15.00 A, a = 89.9°, /? = 109.6°, 
7 = 109.3° and t>o = abc sin a sin f3 sin 7 = 1971 A 3 (see, 
e.g., Ref. 01) O ne can write /3y of Eq. i|12[l in the form 



ft, = 



4£ 



7rA 



V/2 



4(g/i / 5) £ V 2 . (23) 



For Fc 8 £b = 126.4 mK and A ~ 10~ 7 K, so that 

4Ed/(ttA) ~ 1.6 x 10 6 and for not too fast sweep, 

e ~ 10~ 2 , one has £ ~ 10 5 . This is also an estimation 

1/3 

for the number of spins within the distance r c = («o0 
that strongly interact with a given spin. 

Consider a macroscopically large specimen of ellip- 
soidal form. According to Eq. ijTBl) ly does not diverge 
for - — > 00, and for £ ^> 1 one can replace the sum in 



In general, F(x,8 a ) is computed numerically from Eq. 
©■ For 8 a < 1 we use Eq. that yields 

G S -5.73432 + 16 (tt/3) 5/2 |J q | + (l/3 - n (z) ) 4tt. (26) 

This result is non-analytical in 8 a because random inho- 
mogeneities change the asymptotic behavior of F(f3,8 a ) 
at (3 — > —00. The large numerical factor in Eq. I)26|l makes 

G (Sphere) 

very sensitive to 8 a . For 8 a ^> 1 the contribu- 
tion of the integral term in Eq. (|25|l becomes relatively 
small, and one obtains from Eq. (|21|l 



G^ 



K+ (l/3-n( z ))47T 

2%/7T(5 a 



(27) 



Eq. by an integral converging at 7"j 



that makes 



Whereas G < for the sphere and disc, DDI acting 
predominantly antiferromagnetically and enhancing LZ 
transitions, the result for the needle in Eq. H24|) be- 
comes positive already for 8 a > 0.1 (see Fig. P). As 
«i - (VVA)£ 1/2 [see Eqs. JT^ and JSJ], already reso- 
nance shifts Vt of order A that can stem from different 
sources yield a, ~ 8 a ~ 0.1 for sweep rates e ~ 10~ 2 . 
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Let us compare our results with MFA result for S a = 0. 



r(MFA) _ 



'0 



D, 



D zz £, D zz = J2 4>ij 

3 

^(Sphere) + ^1/3 _ „(*)) 47,-. (28) 



Unlike /q in the limit £ ^> 1, the value of /, 



(MFA) 




de- 



pends on the lattice structure. For a simple cubic lat- 
tice D^ phcrc - ) = and the result for D zz becomes purely 
macroscopic. For tetragonal lattices D z S z phcrc ^ > if 
a = b > c and D zz pheTC ^ < if a = b < c. Direct nu- 
merical calculation yields £)(| phcrc ' — 5.139 for Mni2 and 
4.072 for Fe 8 . Note that E = -(1/2)D ZZ E D is the dipo- 
lar energy per site for the ferromagnetic spin alignment. 
Our result E = — A. 131 Ed for the needle-shaped Fegis 
in qualitative accord with Eq = — 4.10Ed of Ref. [12|. 
One can see that /^ MFA ) > q f or the needle and sphere 



whereas L 



(MFA) 




< for the disc, in contradiction with 



our rigorous results above. 

We have shown that the DDI generate huge corrections 
to the standard LZ picture, Iq ~ £ 3> 1, because of its 
long-ranged character. Strong nearest-neighbor interac- 
tion generate only moderate values of Iq, e.g., Io = z/2 
for the ferromagnetic coupling [see Eq. (|14|) ]. and they 
cannot compete with the DDI. Our main result, Eq. J7|), 
is applicable for e|Jo| < 1 so that the term e 2 Io is a 
correction to the leading term e that defines the small 
transition probability 1 — P. The theory breaks down 
in the slow-sweep range e|io| ^ lj where the interac- 
tion strongly modifies the process. Nevertheless Eqs. J7J|, 
(|23[1 . and (|24|) allow to estimate the range of sweep rates 
where the single-particle description of the LZ effect is 
valid and to see whether the interaction tends to suppress 
or to enhance transitions. For the DDI the standard LZ 
effect can be observed for 



(29) 



that for Fe§ results in a rather fast sweep rate e c ~ 2.3 x 
1CP 5 for a sphere without inhomogeneities (G ~ —5.73). 
This would preclude observation of a standard LZ effect 
in experiments. However we have seen above that inho- 
mogeneities of individual resonances drastically reduce 
the effect of interaction and thus increase s c . 

In the sweeping experiments Q on the ±S transi- 
tions in Fes the sweep rate was v = 2Sgfj, B dB /dt, so 
that with Eq. (PJ) one obtains v c = nA 2 /(2he c ) and 
{dB/dt) c = nA 2 /(AHSgu B £ c ) ~ 8 x !Q- r °/e c ~ 3 T/s. 



However Fig. 2 of Ref. |3j shows that (i) standard LZ 
effect can even be seen down to (dB/dt)° c xp ~ 0.01 T/s 
(i.e., e° xp ~ 10 -2 ) and (ii) that A is underestimated for 
dB/dt < (dB/dt) c c xp . This suggest that transitions are 
suppressed, i.e., ferromagnetic couplings are dominating 
and Iq > 0. As in Ref. [3j a crystal of rectangular shape 



(l a = 80 /im, If, — 50 jLtrxi, l c = 10 13:| ) was used 
whose shape is closer to the needle than to the sphere, 
the sign of the effect could be reconciled with our theory 
by assuming even small random inhomogeneities. On the 
other hand, for this sample the inhomogeneities Vi of Eq. 
(|12|l are of the order of the dipolar field itself that leads 
to much larger values of e c , i.e., to slower sweep rates at 
which the interaction precludes observing the standard 
LZ effect. These inhomogeneities are not random, how- 
ever, and the resonance shifts 7^ increase with the dis- 
tance between i and j, depending on the gradient of the 
dipolar field. Dealing with this case would require using 
Eq. J7J) with a complete solution for the inhomogeneous 
dipolar field in the sample. While it can be done else- 
where, we recommend to perform experiments on crys- 
tals of elliptic shape to avoid strong inhomogeneities that 
suppress the effect of interaction and to see a more dra- 
matic influence of the DDI on the LZ transitions. 

To conclude, quantum transitions in a system of many 
interacting two-level particles is a tough problem that in 
general does not yield to familiar approximate methods 
such as the mean-field approximation. We have calcu- 
lated rigorously the staying probability P of the LZ effect 
in the fast-sweep limit, e <C 1, for general interactions. 
We have shown that long-range interactions exceeding 
the level splitting A, such as the DDI, exert a profound 
influence upon the process. For spherical samples DDI 
acts antiferromagnetically (contrary to the MFA predic- 
tions for Mni2 and Fes) and enhances transitions. This 
should lead to overestimating of the experimental value 
of A, if the standard LZ formula, Eq. @, is used. To the 
contrary, for s > 1 one can expect P > e~ £ , determined 
by those particles that are prevented from tunneling by 
positive couplings. 
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